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We study the response of a magnetic-field-driven superconducting qubit strongly coupled to 
a superconducting coplanar waveguide resonator. We observed a strong amplification/damping 
of a probing signal at different resonance points corresponding to a one and two-photon emis¬ 
sion/absorption. The sign of the detuning between the qubit frequency and the probe determines 
whether amplification or damping is observed. The larger blue detuned driving leads to two-photon 
lasing while the larger red detuning cools the resonator. Our experimental results are in good 
agreement with the theoretical model of qubit lasing and cooling at the Rabi frequency. 


Motivated by the first experiment demonstrating the 
energy exchange between a strongly driven superconduct¬ 
ing qubit and a resonator at the Rabi frequency OrD, 
Hauss et alP elaborated a theoretical model to quantify 
this phenomenon. Their model predicts large resonant 
effects for the one- and two-photon resonance conditions 
Qr = uj r — gUh and ULr = 2c o r — g^n, where uj r is the fun¬ 
damental frequency of the resonator, g% is the effective 
coupling energy, and n is the average number of photons 
in the resonator at frequency uj r . Depending on the de¬ 
tuning between the driving frequency UJd and the qubit 
eigenfrequency uj q , either a lasing behavior (blue detun¬ 
ing UJd - ujq > 0 ) of the oscillator can be realized or the 
qubit can cool the oscillator (red detuning UJd - u q < 0 ). 
According to the theory, one-photon lasing/cooling ef¬ 
fects vanish at the symmetry (degeneracy) point of the 
qubit. However, the two-photon processes persist at the 
symmetry point where the qubit-oscillator coupling is 
quadratic and decoherence effects are minimized. There, 
the system realizes a single-atom-two-photon laser. Note 
a similar two-photon lasing by a quantum dot in a mi¬ 
crocavity, which was investigated theoretically in Ref. [3} 

Experimentally, a single qubit one-photon lasing was 
demonstrated by the NEC group 4 . Here, a single charge 
qubit was used and a population inversion was provided 
by single-electron tunneling. Later on, the amplifica¬ 
tion/deamplification of a transmitted signal trough a 
coplanar waveguide resonator was achieved by a strongly 
driven single flux qubitP However, the two-photon lasing 
has not been experimentally demonstrated yet. In this 
paper, we demonstrate the two-photon lasing, as well as 
considerable enhancement of one-photon lasing of a su¬ 
perconducting qubit by one order of magnitude in com¬ 
parison with Ref. [5] This enhancement was achieved by 
a much stronger coupling of the superconducting qubit 
to the resonator. 

The lasing effect was investigated by making use of a 


standard arrangement: a superconducting qubit placed 
in the middle of a niobium A /2 coplanar waveguide 
resonator. The latter was fabricated by conventional 
sputtering and dry etching of a 150-nm-thick niobium 
film. The patterning uses an electron beam lithography 
and a CF 4 ion etching process. The aluminum qubits 
were fabricated by the shadow evaporation technique. 
The coupling between the qubit and the resonator was 
implemented by a shared Josephson junction (Fig. [I]). 
The dimensions of the qubit’s Josephson junctions are 
0.2 x 0.3 /im 2 , 0.2 x 0.2 /im 2 and 0.2 x 0.3 /im 2 , the crit¬ 
ical current density is about 200 A/cm 2 , and the area 
of the qubit loop is 5 x 4.5 /im 2 . The resonance fre¬ 
quency and the quality factor of the resonator’s funda¬ 
mental mode taken for a weak probing (-141 dBm) are 
uj r = 2tt x 2.482 GHz, Qo = 18 000. The same parame¬ 
ters of the third harmonics taken at the same power are 
uj r s = 27 r x 7.446 GHz, Q 3 = 3750. These values were 
determined from the transmission spectra of the coplanar 
waveguide resonator. 
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FIG. 1. a) Scanning electron microscope image of the qubits 
incorporated into the coplanar waveguide resonator, b) De¬ 
tailed scheme of the qubits. The qubits share a Josephson 
junction with each other as well as with the resonator. 


In practice, we measured a two-qubit sample which 
represents a unit cell of a one-dimensional array of fer- 
romagneticalfy coupled qubits exhibiting a large Kerr 
nonlinearity . 6 However, by applying a certain energy bias, 
one qubit can be set to a localized state, while the second 
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is in the vicinity of its degeneracy point. This way, we 
can measure the qubits separately to reconstruct their 
parameters 7 , and the dynamics of the system is defined 
by a single qubit only. Therefore, to describe our find¬ 
ings, we will use the one-qubit model elaborated in Ref. 0 
in which the corresponding Hamiltonian reads in the flux 
basis of the qubit: 

H = - A a x - hCl R0 cos (u> d t)a z (1) 

+ hoj r a^a + gcr z (a) + a) 

where A is the energy level separation of the two level 
system at zero energy bias e = 0 , Hro is the driving 
amplitude of the applied microwave magnetic flux with 
frequency uJd, and g is the coupling energy between the 
qubit and the resonator. The coupling energy scales with 
the ratio of the magnitude of the persistent current in the 
qubit I q and the critical current of the coupling Joseph- 
son junction I c o as 


_ _hu r I q I I 
9 ~ 9 ~ 2 tt J c0 V G 0 Z r 


( 2 ) 


where Z r = 50f2 is the wave impedance of the coplanar 
waveguide resonator and Go = 2 e 2 /h is the quantum 
conductance. 

This Hamiltonian can be transformed by the Schrieffer- 
Wolff transformation U = exp (iS) with the generator 
S = ( g/hujq ) cosr](a + o))<J y and a rotating wave approx¬ 
imation Ur = exp(— iu)d(J z t/2) to the Hamiltonian 2 

H = hu: r a) a + -hI}R<j z 

+ g sin r] [sin / 3cr z — cos /3a x \ (a + aJ) 

g 2 

— ——cos 2 r)[sm/3(r z — cos/3cr x \(a + a ^) 2 

hidq 

Here, hu q = \J e 2 + A 2 , Hr = cos2 V + &^ 2 > 

tan g = e/A, tan/3 = 5lj/(Qro cos 77 ) and Suj = ujd — uj q . 
The transmission of the resonator 


t oc (a) (3) 

where 


(a) = tr(pa) 


(4) 


was calculated numerically by the quantum toolbox 
QutipP solving the Liouville equation for the density ma¬ 
trix of the system in the rotating frame 


P=- 


h - 


H,p 


+ Lqp + L r p 


(5) 


where L q and L r are Lindblad superoperators 

f 

L q p = y ( 2c7 -/ 5cr + - P°+°- ~ °+°-P) 

f t 

+ y ( 2 a+pa- - p<J-a+ - <J-<J+p) 

+ y {ozpoz-p), (6) 


L r p = Nth + l)( 2 apa^ — pa) a — a) ap) 

+ ^N t h( 2 a)jkb - aa)p - paa f ). ( 7 ) 

Here Nth = 1/ [exp(hujT/ksT) — 1] is the thermal distri¬ 
bution function of photons in the resonator, r is resonator 
loss rate, T^ and T^ are the relaxation, excitation and 
dephasing rates in the rotating frame derived in Ref. 2 


ffd = 7^ cos2 vO- i ski P ) 2 + sin 2 77 cos 2 /3 

F 

f ^ cos 2 77 cos 2 P + T^ sin 2 77 sin 2 ft (8) 



FIG. 2. (a) Resonator transmission as a function of the 

magnetic flux and the detuning of the resonator from the third 
harmonic resonance frequency cj r 3/27r. The anticrossings of 
qubit A and B are separated by a magnetic flux of about 
55 x 10 _3 $o- (b) Close view of the transmission in the vicinity 
of the left qubit’s (A) degeneracy point and (c) the cut of the 
transmission map along the dashed line in (b). The crosses 
are experimental data and the solid line is a theoretical curve 
calculated from Eq. [ 5 ] 

The qubit parameters used for the numerical calcula¬ 
tions were determined independently from the transmis¬ 
sion of the resonator t coupled to the undriven qubit. For 
a weak microwave signal with frequency c j s , the transmis¬ 
sion can be expressed in simple form 9 


^ ^ l^ext SuJq + iq /g\ 

2 uo 2 cos 2 77 — (Scj r + ik/ 2) (Sujq + ij) 

where 6cj q = cu q — cj s , Scu r = uj r — uj s , R ext is the ex¬ 
ternal loss rate of the resonator and 7 is the qubit deco¬ 
herence rate. The experimental data was fitted by Eq! 
(see Fig. [ 2 ]) and the qubit parameters obtained from the 
fitting procedure are given in the table [I| 

We have investigated the stimulated emission effect ob¬ 
served when strongly driving the system at a frequency 
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Qubit I p A/27T g/27r 7/27r 

(nA) (GHz) (MHz) (MHz) 

A 208 6.39 109 15 

B 138 5.28 77 20 

TABLE I. Qubit parameters determined from the fitting pro¬ 
cedure. 
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FIG. 3. a) Resonator transmission in dB units as a function 
of the magnetic flux and the detuning of the resonator from 
the resonance frequency 2.482 GHz. 


cJd/27r = 9c<; r /27r = 22.338 GHz for qubit A. The res¬ 
onator transmission was measured by a network analyzer 
at resonance uo r for magnetic fluxes marked by the black 
rectangular area in Fig. [ 3 ] and is shown in Fig. [ 4 ^i. At 
a driving power = —103 dBm, two emission peaks 
(e.l, e.2) accompanied by two attenuation dips (a.l, a.2) 
appear in the transmission spectra. The increase of the 
transmission is accompanied by a narrowing of the res¬ 
onance curve. The one-photon (a.l, e.l) and two pho¬ 
ton (a.2, e.2) processes are enhanced at resonance with 
the Rabi frequency of the qubit Q# = uo r — g%n and 
Q r = 2c u r — g^fi, respectively. These results are in good 
agreement with the theoretical model 2 described above 
for parameters given in Table |IJ By a strong coupling 
of the qubit to the resonator, we have achieved a con¬ 
siderable enhancement of the lasing, nearly one order of 
magnitude, in comparison with the results presented in 
Ref. 0 Further improvement is possible by increasing 
the relaxation rate of the qubit, for instance, by placing 
a gold resistor close to the qubit loop. 

To conclude, we have experimentally demonstrated 
single-qubit one-photon and two-photon lasing. The ex¬ 
perimental results are in good agreement with the theo¬ 
retical model developed by Hauss et alP The considerable 
enhancement of lasing effect was achieved by stronger 
coupling of the superconducting qubit to the resonator, 
and theoretical calculations show that it can be enhanced 
further by increasing the relaxation rate of the qubit. 
Such improvement could enable to observe even higher- 
order processes analysed in Ref. HOl 
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FIG. 4. Transmission of the resonator at fixed frequency 
LJs/Z'K for a driving signal with frequency cjd/ 2 tt = 9o; r /27r = 
22.338 GHz and power pa = —102 dBm for the driving 
switched off and on (a). Regions e.l, e.2 and a.l, a.2 ex¬ 
hibit amplification and attenuation of the signal, respectively. 
Panel (b) shows the simulated average photon number in the 
resonator obtained for the qubit parameters. 
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FIG. 5. Resonance curve of the resonator at driving switched 
off (solid line) and on (dashed line) with frequency uoa/2 tt = 
9o; r /27r = 22.338 GHz and power pd — —101 dBm in region 
e.l. The amplitude increases by a factor of ~ 9 and band¬ 
width is reduced by factor of 10. 
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